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Abstract. We investigate the axisymmetric structure of col- 
limating, relativistic, strongly magnetized (force-free) jets. In 
particular, we include the differential rotation of the foot points 
of the field Unes in our treatment. The magnetic flux distribu- 
tion is determined by the solution of the Grad-Shafranov equa- 
tion and the regularity condition along the light surface. With 
differential rotation, i.e. the variation of the iso-rotation param- 
eter rip, the shape of the light surface is not known a priori and 
must be calculated in an iterative way. 

For the first time, we have calculated the force-free mag- 
netic structure of truly two-dimensional, relativistic jets, an- 
chored in a differentially rotating disk. Such an approach al- 
lows for a direct connection between parameters of the central 
source (mass, rotation) and the extension of the radio jet. In 
particular, this can provide a direct scaling of the location of 
the asymptotic jet light cylinder in terms of the central mass 
and the accretion disk magnetic flux distribution. 

We demonstrate that differentially rotating jets must be col- 
limated to a smaller radius in terms of the light cylinder if 
compared to jets with rigid rotation. Also, the opening angle 
is smaller Further we present an analytical estimate for the jet 
opening angle along the asymptotic branches of the light sur- 
face. In general, differential rotation of the iso-rotation param- 
eter leads to an increase of the jet opening angle. 

Our results are applicable for highly magnetized, highly 
collimated, relativistic jets from active galactic nuclei and 
Galactic superluminal jet sources. Comparison to the M87 jet 
shows agreement in the coUimation distance. We derive a light 
cylinder radius of the M87 jet of 50 Schwarzschild radii. 

Key words: Accretion, accretion disks - MHD - Methods: 
numerical - ISM: jets and outflows - Galaxies: individual: M87 
- Galaxies: jets 



1. Formation of magnetic jets 

Observations of astrophysical jet sources have now established 
the idea that jet formation is always connected to the pres- 
ence of an accretion disk and strong magnetic fields. This holds 
for various scales of energy output, jet velocity and nature of 
the jet emitting objects. Examples are jets from active galac- 



tic nuclei (AGN), Galactic superluminal jet sources, the ex- 
ample of a mildly relativistic jet from a neutron star (SS 433) 
and the numerous class of protostellar jets (see Zensus et al. 
1995; Mirabel & Rodriguez 1995; Mundt et al. 1990; Ray et al. 
1996). Magnetic jets are believed to originate very close to the 
central object in the interaction region with the accretion disk. 
Beside observational arguments also theoretical considerations 
have shown that magnetic fields play an important role in jet 
formation and propagation (Blandford & Payne 1982; Pudritz 
& Norman 1983; Shibata & Uchida 1985; Sakurai 1985; Ca- 
menzind 1987; Lovelace et al. 1991). 

If the central object is a black hole as it is the case for AGN 
and Galactic superluminal jet sources, the surrounding accre- 
tion disk is the only possible location for a field generation (by 
dynamo action or/and advection of flux). In the case of stellar 
objects (protostars, white dwarfs or neutron stars), the central 
star also carries a relatively strong magnetic field and it is not 
yet clear, whether the jet magnetosphere originates in the disk 
or in the star. However, a strong interaction between stellar field 
and accretion flow is evident. The jet formation process itself is 
not yet fully understood theoretically. In particular, for the mass 
transfer from the disk into the jet and the process of magnetic 
field generation a complete physical model is missing. 

However, over the last decades the basic ideas of Blandford 
& Payne (1982) have been extended by various authors. The 
general picture is the following. Matter is lifted from the disk 
into the magnetosphere and becomes magnetically accelerated 
(Ferreira 1997). Toroidal magnetic fields, generated by inertial 
back-reaction of the plasma on the poloidal field, may colli- 
mate the disk magnetosphere into a highly collimated jet flow 
(Camenzind 1987;Chiuehetal. 1991; Lovelace et al. 1991). In 
general, due to the complexity of the MHD equations, station- 
ary relativistic models of magnetic jets has to rely on simpli- 
fying assumptions such as self-similarity (Contopoulos 1994, 
1995), some other prescription of the field structure (Li 1993, 
Beskin 1997) or the restriction to asymptotic regimes (Chiueh 
et al. 1991; Appl & Camenzind 1993; Nitta 1994, 1995). For 
highly magnetized jets the force-free limit applies. This allows 
for a truly two-dimensional calculation of the magnetic field 
structure (Fendt et al. 1995; Fendt 1997a). Relativistic magne- 
tohydrodynamics implies that poloidal electric fields, which are 
not present in Newtonian MHD, cannot be neglected anymore. 
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From the observations we know that extragalactic jets as 
well as Galactic superluminal jets and protostellar jets are col- 
limated almost to a cylindrical shape (Zensus et al. 1995; Ray 
et al. 1996; Mirabel & Rodriguez 1995). Theoretically, it has 
been shown that cuiTent carrying relativistic jets must collimate 
to a cylinder (Chiueh et al. 1991). For the asymptotic limit of 
a cylindrically collimated magnetic relativistic jet, Appl & Ca- 
menzind (1993) found a non-linear analytical solution for the 
trans-field force-balance in the case of a constant iso-rotation 
parameter These results were further generalized for jets with 
differential rotation (Fendt 1997b). Such an asymptotic field 
distribution is especially interesting for jets emerging from (dif- 
ferentially rotating) accretion disks. 

In previous papers, we applied the asymptotic jet model 
of Appl & Camenzind (1993) as a boundary condition for the 
calculation of gZo^aZ, two-dimensional, force-free jet magneto- 
spheres for rapidly rotating stars (Fendt et al. 1995) or rapidly 
rotating black holes (Fendt 1997a). In this paper, we continue 
our work on 2D force-free jet magnetospheres applying an 
asymptotic jet with differential rotation of the iso-rotation pa- 
rameter flp as boundary condition for the global jet structure. 
Such an approach allows for a connection between the differen- 
tially rotating jet source - the accretion disk - and the asymp- 
totic collimated jet. Since jet motion seems intrinsically con- 
nected to the accretion disk, differential rotation of the field 
lines should be a natural ingredient for any magnetic jet struc- 
ture. As a principal problem for differentially rotating relativis- 
tic jet magnetospheres, position and shape of the singular light 
surface are not known a priori, but have to be calculated in an 
iterative way together with the magnetic flux distribution. 

In Sect. 2 we recall some basic equations of the theory of 
relativistic magnetospheres and discuss several difficulties with 
the solution of the Grad-Shafranov (hereafter GS) equation. 
After some comments on the numerical approach in Sect. 3, we 
discuss our results in Sect. 4. A summary is given in Sect. 5. 



which, in the force-free case, flows parallel to the flux surfaces, 

The structure of the magnetic flux surfaces is determined by 
the toroidal component of Ampere's law, V x Bp = Airj^/c, 
where the toroidal electric current density has to be calculated 
from the equation of motion projected perpendicular to the flux 
surfaces (Camenzind 1987; Fendt et al. 1995). For strong mag- 
netic fields, inertial forces of the matter can be neglected. This 
is the force-free limit and the equation of motion reduces to 
= cpcE + j X B with the charge density p^. 

Combining Ampere's law and the force-free equation of 
motion the cross-field force-balance can be written as the mod- 
ified relativistic GS equation. 



4 1 1 
'^R2 



(/'(*))' (3) 



where the primes denote the derivative ^ (see Camenzind 
1987; Okamoto 1992). 

fip is conserved along the flux surfaces, ftp = f7F(^')- We 
will call it the iso-rotation parameter, defined by Ferraro's law 
of iso-rotation. It can be understood as the angular velocity of 
the plasma, reduced by the slide along the field lines. Some- 
times, it is called the angular velocity of the field lines. Both, 
the current distribution /(\E') and the rotation law Of(\E') de- 
termine the source term for the GS equation and govern the 
structure of the magnetosphere. We apply the following nor- 
malization, 

R, Z xRq,zRq, 

rip ^ f^F (c/i?o) , 
I ^ I /max ■ 



2. Structure of magnetic jets 

Throughout the paper we apply the following basic assump- 
tions: axisymmetry, stationarity and ideal MHD. We use cylin- 
drical coordinates (i?, 0, Z) or (x, (/), z) if normalized. The 
term 'asymptotic' always denotes the limit of Z » R un- 
less explicitly stated otherwise. We consider jets with a. finite 
radius, R < oo for Z oo. 



2.1. The force-free, cross-field force-balance 

With the assumption of axisymmetry, a magnetic flux function 
can be defined 



* = -5- / Bp • dA, i?JBp = V* A 



(1) 



measuring the magnetic flux through a surface element with 
radius R and, in a similar way, the poloidal current through the 
same area 



jp • dA = --RB^ 



(2) 



As the length scale for the GS equation (3) the asymptotic ra- 
dius Rq of the light surface is selected (see below). In order 
to allow for a direct comparison to rigidly rotating magneto- 
spheres, the normalization was chosen such that Jlp = 1 at 
X = 1. With the chosen normalization, Eq. (3) can be written 
dimensionless. 



i-x'^nK^) 



1.9 



- x\w^\^^{nU^)y. (4) 

5 is a coupling constant describing the strength of the current 
term in the GS equation. 



^ _ ^-'max^O _ ^ / ^max \ / J^O 



^ ^ max 
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where the parameters are chosen for extragalactic jets. In this 
paper, g is in accordance with the definition in Fendt et al. 
(1995) and differs from the definition in Appl & Camenzind 
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(1993) by a factor of two, gpcndt = 2 (^acQ Interestingly, a 
coupling constant, defined in a similar way also for the differ- 
ential rotation term, would be equal to unity. The GS equation 
is numerically solved applying the method of finite elements 
(see Appendix). 

Along the Ught surface, where D = 1~ x^fl'^{'$) — 0, the 
GS equation reduces to the regularity condition, 

V* . Vi^ = -g^ - ilV^/p (In (17f(*)2))' , (5) 

which is equivalent to a Neumann boundary condition. How- 
ever, for differentially rotating magnetospheres with Hp = 
r2F(^') the shape of this surface is not known a priori and 
has to be calculated in an iterative way together with the two- 
dimensional solution of the GS equation. For constant Op the 
light surface is of cylindrical shape. As we have shown in a pre- 
vious publication (Fendt et al. 1995), our finite element code 
satisfies the regularity condition automatically, since the sur- 
face integral along the light surface vanishes. 

2.2. Discussion of the force-free assumption 

It is clear that relativistic jets must be highly magnetized. Only 
a high plasma magnetization gives jet velocities close to the 
speed of light (Fendt Sl Camenzind 1996). Therefore, for the 
calculation of field structure the force-free limit seems to be 
reasonable. However, one may question the assumption of a 
force-free asymptotic jet. In a fully self-consistent picture of 
magnetic jet formation, the asymptotic jet is located beyond 
the coUimating, non force-free wind region and beyond the fast 
magnetosonic surface. The asymptotic jet parameters are de- 
termined by the wind motion. Thus, poloidal current and iso- 
rotation parameter of the field are not functions free of choice. 
The force-free region of a jet is located close to its origin, where 
the speed is low. Beyond the Alfven surface plasma kinetic en- 
ergy dominates the magnetic energy, which is just the contrary 
to force-freeness. 

For small plasma density, the Alfven surface of the wind 
flow approaches the light surface. In this case the fast magne- 
tosonic surface moves to infinity for a conical flow. Okamoto 
(1999) argues that a force-free field distribution extending to 
infinity in both x and z direction will asymptotically be of con- 
ical shape, i.e. un-collimated. However, his approach differs 
from ours in the sense that he assumes that all field lines will 
cross the light cylinder Such an assumption per se prohibits 
any collimation. On the other hand, perfect jet collimation is 
an observational fact. Astrophysical jets (of very different en- 
ergy scales) appear collimated to cylinders of finite radius. 



' Due to the fact that the jet radius (where ^' = 1) is not Imown 
before the asymptotic GS equation has been solved (Fendt 1997b), the 
normalization with g leads to a current distribution which is not 
normalized to unity. This difference in normalization is "hidden" in 
the coupling constant g, which could, in principal, be re-scaled appro- 
priately. 



In general, the non force-free relativistic GS equation 
shows three inertial contributions, 

= -i (1 - M-' - I'nl) + {1 - x^nf) ^ + ^ 

where Vl indicates the gradient perpendicular to k = 
kBp/Att = n • (Bp • V)Bp/47r the poloidal field curvature, 
p the mass density, the toroidal velocity, P the gas pressure 
and M the Alfven Mach number (Chiueh et al. 1991). One can 
show that in the asymptotic, cylindrical jet the contribution of 
inertial terms in the force-balance across the field is weak. The 
contribution of gas pressure is usually negligible in astrophysi- 
cal jets. Also, the centrifugal force does not play a role for radii 
larger than the Alfven radius, since outside the Alfven surface 
(where = 1 — x^ilp) the plasma moves with constant an- 
gular momentum. The curvature term vanishes in cylindrical 
geometry. Therefore, since for cylindrical jets the contribution 
from inertial terms is weak, the configuration is comparable to 
the force-free case. The source term of the GS equation may be 
reduced to a form similar to the common force-free limit. We 
suggest the phrase "quasi force-free" for such a configuration 
because the GS equation looks force-free even if the physical 
system is not magnetically dominated. 

In the force-free limit of a highly magnetized plasma the 
previous arguments also apply. However, in difference to the 
asymptotic regime considered above, the low plasma den- 
sity implies that inertial terms are weak over the whole two- 
dimensional jet region. The centrifugal term pu^ is weak even 
if the Alfven surface now comes close to the light surface. Nu- 
merical calculations of the plasma motion along the field have 
shown that, for a high magnetization, the Alfven Mach num- 
ber M grows almost linearly with radius but remains relatively 
low (Fendt & Camenzind 1996). Thus, the inertial curvature 
term should not play a dominant role. Contopoulos & Lovelace 
(1994) find from self-similar solutions that centrifugal forces 
are dominated by magnetic forces leading to a re-collimation 
of the outflow. 

In summary, our discussion of the inertial terms in the 
force-balance equation has shown that these terms are gener- 
ally weak in the case of a high magnetization. We therefore 
think that for the calculation of the magnetic field structure 
in relativistic jets the force-free assumption is acceptable. The 
main motivation of the force-free assumption is clearly the rea- 
son of simplification. There is yet no other way to calculate 
a truly two-dimensional field distribution for relativistic jets. 
Naturally, with a force-free solution, nothing can be said about 
the flow acceleration itself. 

2.3. Location of the asymptotic light cylinder 

The radius of the asymptotic light cylinder Rq is the natural 
length-scale for the GS solution. A scaling of the GS solution 
in terms of the properties of the central object (e.g. its mass) 
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relies on the proper connection between the asymptotic jet and 
the accretion disk. This is feasible only if differential rotation 
riF(^') is included in the treatment (see Sect. 3). 

In the following we consider the location of the light sur- 
face and its relation to the relativistic character of the mag- 
netosphere from a general point of view. The light cylinder 
of a flux surface 5* is defined as a cylinder with radius R = 
i?L(^) = c/^^f(^)- At this position the GS equation becomes 
singular However, this light cylinder is only important for the 
field line ;/ the field line actually intersects it as for '^out in 
Fig. 1. Only then, relativistic effects become dominant. For ex- 
ample, the poloidal electric field scales with the radius in units 
of the light cylinder radius, E-p = {R/ R\,)Bp. On the other 
hand, in the case of 'fin in Fig. 1, the asymptotic radius of the 
flux surface is smaller than its light cylinder radius i?L(^in) 
(located between ^'in and 'I'o), therefore relativistic effects are 
small. For jet solutions with rigid rotation Vl-p all flux surfaces 
have the same light cylinder radius. Thus, the singular light sur- 
face of the magnetosphere is a cylinder For jet solutions with 
differential rotation ilp the flux surfaces have different light 
cylinder radii. The singular surface of the magnetosphere is not 
a cylinder anymore. 

It is now interesting to note that the case of differential ro- 
tation fiF(5') allows for a hypothetical field distribution where 
(i) the light radius of most of the flux surfaces is located within 
the jet radius, but where also (ii) the asymptotic radius of the 
flux surfaces is always smaller than their light radius. Such a 
field distribution would not have a singular light surface and 
could be considered as "non relativistic", even if the hypotheti- 
cal light radii of many field lines are inside the jet radius. Such a 
situation is not possible if the magnetosphere is constrained by 
a constant rotation Vl-p- This underlines the importance of the 
treatment of differential rotation for jets from accretion disks. 

A relativistic description for the jet magnetosphere is al- 
ways required if the jet contains a flux surface for which the 
light radius is smaller than the asymptotic radius. 



asymptotic collimated jet 




Fig. 1. Sketch of the jet model. Axisymmetric jet magnetic flux 
surfaces 5* projected into the meridional plane. The central 
object, located within the inner boundary (solid disk), is sur- 
rounded by an accretion disk. Helical magnetic field lines (lay- 
ing on the flux surfaces) are anchored in the differentially rotat- 
ing disk at the foot points The jet boundary is defined 
by the flux surface = 1. The upper boundary condition is 
a cylindrically collimated jet solution (Fendt 1997b). The ar- 
row indicates the numerical deformation of the initially verti- 
cal boundary of the inner solution (at x ^ V) into the curved 
light surface. The flux surfaces ^E'in (4'out) have an asymptotic 
radius smaller (larger) than the asymptotic light cylinder i?o, 
which is the asymptotic branch of the light surface i?L(^') for 
large z. The flux surface vJ/q coincides with the light surface 
asymptotically. The jet half opening angle is a (see Sect. 2.4, 
Fig. 2). 



2.4. The regularity condition and the jet opening angle 

The regularity condition (5) is the natural boundary condition 
along the light surface. Although it is impossible to solve equa- 
tion (5) explicitly, a general relation concerning the jet opening 
angle can be derived. First, we rewrite Eq. (5) as 

B.^i,a')'-iBf(^)', (6, 

where ilpl*) = = R^/Ri^i'^) has been applied. 

From Eq.(6) it follows for the radial field component = 
-g(/^(\E'))7(l/^^|(*))', if * intersects the light surface with 
vanishing ■ On the other hand, considering a field line per- 
pendicular to the light surface, V^* _L VD, this provides a 
condition for the axial field component. 



Interestingly, this is either only a function of the current dis- 
tribution /(^) or depends only from the specification of the 
rotation law Of(^'). Further, in this case it is always B, > 0, 
since (l/fi|)' = (s^l)' > 0. In particular, for the asymptotic 
(z — > oo) part of the magnetosphere, this implies that only col- 
limating field lines can cross the light surface. 

Now we consider the asymptotic branches of the light 
surface. For the asymptotic branch in z-direction it holds 
{VD):, » (V£»)^ ~ 0. Further, it is B^(lnl7|)' = 0, im- 
plying either a collimated field structure, B^; = or rigid rota- 
tion, (f2F(\i/))' = 0. From this we conclude that in the asymp- 
totic regime of a cylindrical light surface, also the flux surfaces 
along this light cylinder must be of cylindrical shape. CoUima- 
tion must occur in the non-asymptotic region of the jet. 

If we now assume that there exists an asymptotic part of 
the light surface in a--direction (where x >> z) and that 
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Fig. 2. Jet half opening angle ol{^ = 1) for the analytical 
asymptotic jet solution and along the asymptotic branch of 
the light surface in x-direction (see Eq. (8)). Coupling constant 
t; = 0.1 (a) and g ^ 2 (b). Asymptotic magnetic flux distribu- 
tion parameter (Eq. (9)) b = 0.5 (solid), 6=1 (dotted), 6 = 2 
(short-dashed), 6 = 5 (long-dashed). 



{yD)z » {S7D)x — 0, we derive an equation for the half 
jet opening angle. 



aW=tan J1 + -.9 



(8) 



for the flux surfaces in this region. As a general example we 
apply the analytical solution obtained for the asymptotic jet 
(Fendt 1997b), 



^{x) = ihi ( 1 



6 = In ( 1 
1 



(l-e-fe*)2 (l-e-'')2 



fjctV 

a / 
+ ^1(1) (9) 



for Eq. (9). Here, 6 is a measure for the ratio of jet radius to jet 
core radius a. Finally, we obtain the half opening angle for the 
outermost flux surface = 1, 



where is defined as {^^{1))' /flp{l). Note, that Eq. (10) is 
only valid in the limit of {WD)x — 0. Fig. 2 shows the vari- 
ation of the opening angle a with the parameters a and 6 for 
two choices of the strength of the poloidal electric current. In 
general, jets with a strong differential rotation 17f(5') (i.e. large 
Rn) have a larger opening angle. Also, jets with a large ratio of 
jet to core radius have a smaller opening angle. Therefore, jets 
originating in a small part of the accretion disk, equivalent to 
small value of Rq, will be collimated to a smaller opening an- 
gle. It is interesting to note that, in the case of a rigid rotation, 
the limiting half opening angle is 45°, independently of g and 
6. 

3. The jet-disk connection, providing the true length scale 
of the GS solution 

The natural length scale of the relativistic GS equation (3) is the 
asymptotic light cylinder Rq (see Sect. 2.3). Its size is related 
to the iso-rotation parameter f2F(^'), which itself is connected 
to the angular rotation of the foot points of the field lines. Con- 
cerning the GS equation, the size of Rq follows purely from 
electro-magnetic quantities, if the coupling constant g is cho- 
sen. The GS solution can be scaled to any central object from 
stars to galactic nuclei as long as the interrelation of the pa- 
rameters ^'max, .^max and i?o providcs the same g. So far, no 
connection has been made to the type of central object. Here, 
we treat the question where the asymptotic light cylinder is lo- 
cated in physical units. 

In the case of rigid rotation, the light cylinder radius is usu- 
ally estimated by choosing a distinct radial distance from the 
central object and defining Vl-p under the assumption that the 
jet magnetosphere is anchored in that point. If the central ob- 
ject is a black hole, the marginally stable orbit impUes an upper 
limit for fip. For jets in AGN this estimate leads to the com- 
mon conclusion that the light cylinder radius is at about 10i?s 
and the typical jet radius at about 100 i?L- Clearly, such argu- 
ments relies on the assumption that the field line emerging at 
this very special radius defining fip also extends to the light 
cyUnder radius i?L (see Sect. 2.3). 

This picture changes, if differential rotation ilF(^) is con- 
sidered. Different flux surfaces anchor at different foot point 
radii and have different Ught radii (Sect. 2.3). Assuming a Ke- 
plerian rotation, the light surface radius i?L(^') is located at 



i?L(^') = 4 X 10^'^ cm 



/ flD(^) 



3/2 



M 



IQio Mr. 







(11) 



where i?g is the Schwarzschild radius of a point mass and 
i?D ('f) the foot point of the flux surface 'I' on a Keplerian disk. 
A more general equation is 



i?L(£) 



= V2 



3/2 



(12) 



tan 



1 + i?n 4i?a 



1 



he'' 



5/2 



i?D ('I') is determined from the magnetic flux distribution along 
^(10) the disk and is defined by a certain disk model. Fig. 3 shows 
the location of the light radius i?L for a field line anchored at 



6 



Fendt & Memola: Relativistic magnetic jets 



100 



K 
03 




2 4 6 

log M (solar masses) 



10 



Fig. 3. Location of the light cyhnder radius of a flux surface 
i?L(4'), anchored at a certain foot point radius i?D(^) in units 
of the Schwarzschild radius i?g in a Keplerian disk around a 
point mass M. Note that for non collapsed stellar objects the 
footpoint radii of the jet field fines are located at about 10*^i?g. 

a foot point radius i?D in a Keplerian disk around a central 
object of mass M (see Eq. 11). Note that the unit of the field 
line footpoint radius in Eq. (12) and Fig. 3 is the Schwarzschild 
radius. Therefore, Fig. 3 is appropriate only for relativistic jets. 
The footpoint radii for protostellar jets are several stellar radii, 
corresponding to about 10^ Schwarzschild radii (which would 
be located far above the box in Fig. 3). 

So far, nothing can be said about the location of the asymp- 
totic radius of the field lines in general. The essential question 
is where the asymptotic radius of a flux surface is located in re- 
spect to its light cylinder This question can only be answered 
by a detailed model considering the two-dimensional field dis- 
tribution including differential rotation Of(^'). Only in such a 
model, the flux distribution of the asymptotic jet can be con- 
nected to the flux distribution of the "star"-disk system. Cer- 
tainly, both boundary conditions - asymptotic jet and accretion 
disk - rely on model assumptions. However, in a self-consistent 
model these boundary conditions have to follow certain con- 
straints (see Sect. 4.1, 4.3). 

4. The two-dimensional jet solution 

4.1. Disk and jet boundary condition 

Three important boundary conditions determine the two- 
dimensional flux distribution. The first boundary condition is 
in the asymptotic region. Here we assume a cylindrically col- 
limated jet. We apply the magnetic flux distribution derived 
by Fendt (1997b), where the rigidly rotating jet model of 
Appl & Camenzind (1993) is extended for differential rota- 
tion 17f(5'). In particular, our asymptotic jet shows the typ- 
ical jet core-envelope structure of magnetic flux and electric 
current, i.e. a configuration where most of the magnetic flux 



and poloidal electric current is concentrated within a "core" 
radius. The asymptotic model provides not only the asymp- 
totic magnetic flux boundary condition but also the flp (^') and 
7(5') distribution for the whole two-dimensional jet magne- 
tosphere. In the model of Fendt (1997b) these functions fol- 
low from the solution of the one-dimensional (asymptotic) 
GS equation across the cylindrical jet and the prescription of 
I{x) ^ (a;/a)V(l + (x/af) together with nl{x) = e''"''^, 
where a is the core radius of the electric current distribution 
and h governs the steepness of the fip-profileQ. 

The second boundary condition is the magnetic flux distri- 
bution along the disk. This distribution is in general not known 
as a solution of the full MHD disk equations. Typical models 
rely on various simplifying assumptions, as e.g. stationarity, the 
distribution of magnetic resistivity or the disk turbulence gov- 
erning a dynamo process. We apply an analytic flux distribu- 
tion similar to the model of Khanna & Camenzind (1992), who 
calculated the stationary accretion disk magnetic field structure 
around a super massive black hole. The typical behavior of the 
magnetic flux distribution is (i) a slight increase of magnetic 
flux along the innermost disk, (ii) a small or vanishing flux at 
the inner disk radius, (iii) a strong increase of magnetic flux 
at intermediate radius (the core radius) and (iv) a saturating 
behavior for large radii. Using the normalization introduced 
above, we choose the following disk boundary magnetic flux 
distribution. 



*disk(2;) = -^In ( 1 




(13) 



with b = \n{l+{xdisk—Xia) / a )) (see Fig. 4). The parameters 
are: the core radius 5, the disk outer radius, I'disk and the disk 
inner radius, Xin- For simplicity, we choose Xin — without 
loss of generality. Such a choice will definitely not influence 
the global jet solution which is normalized to the asymptotic 
light cylinder radius. 

The third boundary condition is the jet boundary Xjct{z)- 
Along this boundary the flux distribution is 5* = 1 by defini- 
tion. However, the shape of this boundary is not known a priori. 
It must be determined by the regularity of the solution across 
the light surface (see also Fendt et al. 1995). A slightly different 
shape may give the same solution. However, the main features 
of the solution as opening angle or locus of the collimation are 
fixed by the internal equilibrium. Therefore, the regularity con- 
dition governs the shape of the jet boundary. For a given 1(5'), 
17f(5'), disk and jet boundary condition, the jet radius Xjct{z) 
is uniquely determined. 

The inner spherical grid boundary with radius x,t close to 
the origin, indicates the regime of the central source, possibly 
enclosing a collapsed object. Neutron stars or magnetic white 
dwarfs may carry their own magnetic field, a black hole can 
only be thread by the disk magnetic field. In any case, the mag- 
netic flux distribution is a combination of "central" magnetic 
flux and disk magnetic flux 4* = + 5'disk- For simplicity we 

" For figures of these functions and the related "^{x), Qf{'^) and 
/(*) distribution, we refer to Fendt (1997b). 
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0.2 0.4 0.6 0.8 



Fig. 4. Magnetic flux distribution along the disk ^'disk(2;) 
as defined in Eq. (13). Parameters; ~ 0, .Tdisk/a = 
100,40,15,7,4,2,1,0.01. 



assume that the magnetic flux increases monotonously from the 
axis to the disk edge and (a;*) ~ 'I'j,(a;^) and 5'disk(a;*) = 0. 

4.2. The two-dimensional collimating magnetic field structure 

Results of numerical solutions of the GS equation are pre- 
sented in Fig. 5. Shown is the two-dimensional structure of the 
magnetic flux surfaces as projection of the helical field lines 
onto the meridional plane. In general, for a choice of the "free 
functions" /('f) and Of(^'), here taken from the asymptotic 
cyUndrical jet solution, the field structure is determined by 
the boundary conditions and the regularity condition along the 
light surface. 

We calculated two solutions with a different choice for the 
steepness parameter in the iso-rotation Hp. The first solution 
is for h = 0.2 (Fig. 5, left). This is more close to the case of 
rigid rotation. Indeed the solution look rather similar to the so- 
lutions presented in Fendt et al. (1995). The second solution 
is for h = 0.5 (Fig. 5, right). The steeper profile for the rota- 
tion law implies a smaller asymptotic jet radius (Fendt 1997b). 
This can be seen in comparison with the rigid rotation solutions 
(Fendt et al. 1995). However, a larger poloidal electric current 
can balance the effect of differential rotation. Therefore, the 
h = 0.2 solution (with g = 2.5) collimates to a smaller asymp- 
totic jet radius than the h = 0.5 solution (with g = 2.0). A 
h = 0.2 solution with g = 2.0 would have a jet radius of 2.4. 
The second solution with the steeper profile of the rotation law 
rjpC^') would better fit to a Keplerian disk rotation. A perfect 
match would require an even steeper riF(5')-profile (see be- 
low). 

The mean half opening angle of the jet magnetospheres is 
about 60°. As discussed above, the shape of the outermost flux 
surface (4* = 1) is not prescribed but is a result of our cal- 
culation eventually determined by the regularity condition. Af- 
ter crossing the light surface the jets collimate to their asymp- 
totic radius within a distance from the source of about 1 — 2 i?o 
along the jet axis. The opening angle of the second solution is 
smaller, however, the jet collimation is achieved only at a larger 



distance from the central source. In our examples, the "jet ex- 
pansion rate", which we define as the ratio of the asymptotic jet 
radius to the foot point jet radius (the "disk radius"), is about 
10. The true scaling of the jet magnetosphere in terms of the 
size of the central object can be determined by connecting the 
jet iso-rotation parameter Of(^') to the disk rotation (see next 
Sect). 

We note that, although in our computations the jet boundary 
Xjct{z) is determined by the force-balance within the jet, and 
therefore subject to the regularity condition, with our results we 
do not prove the magnetohydrodynamic self-collimation pro- 
cess of the jet flow. Clearly, the calculated jet magnetosphere 
is self-collimated in the sense that its structure has been de- 
termined only by the internal properties. However, the actual 
collimation process of the jet flow from an un-coUimated coni- 
cal outflow into a cylinder could only be investigated by time- 
dependent simulations taking into account the interaction with 
the ambient medium. 

On the other hand, we can assume that our jet solution 
is embedded in an ambient medium. If we further assume an 
equilibrium between the internal pressure (magnetically dom- 
inated) and external (gas) pressure along the jet boundary, 
we may derive the gas pressure distribution in the ambient 
medium, since we know the magnetic pressure distribution 
along the collimating jet radius. In this case, the jet solution 
may be considered as collimated by ambient pressure. 

To our understanding one may claim a self-collimation 
only, if the jet flow collimates independently from external 
forces. Since in our treatment we do not consider the interrela- 
tion with the medium outside the jet, we cannot decide whether 
the flow is self-collimated or pressure collimated. 

The field structure is governed by the choice of the func- 
tions /(^') and fipC^'), here taken from an asymptotic jet so- 
lution. In combination with the disk magnetic flux distribution 
(13) we can determine two parameters interesting for the jet- 
disk interaction. These are (i) the magnetic angular momentum 
loss per unit time and unit radius from disk into the jet and 
(ii) the toroidal magnetic field distribution along the disk. With 
as the angular momentum flux per unit time per unit flux 
tube, the (normalized) angular momentum flux per unit time 
per unit radius is dJ/dx = —xB^Iix) along the disk. Fig. 6 
shows the behavior of both quantities for our jet model with the 
steeper profile of the rotation law, = 0.5. As we see, most of 
the magnetic angular momentum is lost in the outer parts of the 
disk. This may have interesting applications for accretion disk 
models taking into account a magnetized wind as a boundary 
condition. The total magnetic angular momentum loss is deter- 
mined by the normalization, ,/ = — (/max* max/ c) / /(^')o?^' 
or j = -(^*max/2i?o) /-/(*)</*■ The magnetic toroidal 
field distribution along the disk has a maximum at about half 
the disk radius. 

Clearly, these parameters are biased by the magnetic flux 
disk boundary condition (13) of our model. However, we be- 
lieve that the main features are rather general and valid for any 
poloidal current and magnetic flux distribution with the typical 
core-envelope structure. 
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Fig. 5. Two-dimensional magnetic flux distribution '^{x, z) for two different asymptotic rotation laws. Left: h ~ 0.2, g — 2.5 
Right: h ~ 0.5, g = 2.0. Shown are iso-contours of the magnetic flux (equivalent to poloidal magnetic field lines) with contour 
levels = lO"*^"'^"-* , n = 0, .., 25. Note that due to the choice of contour levels the iso-contour density does not mirror the 
field strength. 



4.3. Scaling relations of disk and jet 

As discussed above, the two-dimensional magnetic field dis- 
tribution connecting the asymptotic jet region with the lower 
disk boundary allows for a direct scaling of the jet in terms of 
the size of the central object. This is simply based on the as- 
sumption that the foot points of the field lines are rotating with 
Keplerian speed, Vl-p = fix and to the fact that in ideal MHD 
the iso-rotation parameter fip is conserved along the field lines. 
It is therefore possible to construct a self-consistent model of 
the whole "star"-disk-jet system with only a small set of free 
parameters. In the following we will motivate such a model. 

The first example demonstrates how the connection be- 
tween the asymptotic jet and the disk, applied for our very 
special model assumption, provides a specific estimate for the 
asymptotic light cylinder i?o- Normalizing the Keplerian ve- 
locity f^K in the same way as il-p (Sect. 2.1), we obtain the 
expression 

GM 1 ^ GM/c^ ^ 0.5Rs 



Iso-rotation parameter riF(^') and disk radius .Tdisk are fixed 
by our model. Therefore, the asymptotic light cylinder is pro- 
portional to the mass of the central object. For = 0.54 
(which refers to the h = 0.5 model) and .Xdisk = 0.2 the asymp- 
totic light cylinder is Rq — 116i?s, which is about 2 times 
larger compared to the jet solution with a rigid rotation Hp = 1 
and will increase for larger values of h. With the choice of g, 
the value of Rq constraints the maximum poloidal magnetic 
flux and electric current. Here, no assumption is made about 
the flux distribution along the disk. 

In the second example we determine the disk magnetic 
flux distribution 'I'(x) combining the asymptotic jet rota- 
tion law ripl^) with a Keplerian disk rotation Uk{x). From 
Eq. (14) follows that f^Fl^-j/f^Fll) = f7K(a:)/f7F(^' = 1) = 
(a;/a;disk)~'^^^- In combination with the numerically derived 
rip (5') this gives the '^{x) along the disk (Fig. 7). The figure 
shows that the disk flux distribution derived from the asymp- 
totic jet is distributed only over the outer part of the disk. This 
can be interpreted in two ways. First it may imply a relatively 
large inner disk radius and, hence, an asymptotic jet radius 
small in terms of radii of the central object. Secondly, it just 
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Fig. 6. Magnetic angular momentum loss per time unit per unit 
radius dJ/dx at radius x (above) and disk toroidal field distri- 
bution B^{x) (below) for the jet solution with h = 0.5 shown 
in Fig. 5. 

underlines the fact that in our model the distribution of the 
asymptotic jet iso-rotation parameter is too flat in order to be 
truly connected to a disk magnetic flux with an extended radial 
distribution. For a model taking into account the disk Keplerian 
rotation in a fully self-consistent way, the magnetic flux distri- 
bution which has to be used as disk boundary condition for the 
GS solution is the one derived in Fig. 7. 

On the other hand, the assumption of a Keplerian disk rota- 
tion together with a certain disk magnetic flux distribution pro- 
vides an expression for the iso-rotation parameter JIfI^) = 

nKix{^)) = (GA//i?oc2) (e^* - l)) Here, the 
disk magnetic flux distribution (13) has been used. Eventually, 
one finds 



fiF(l) 




-3/4 



(15) 



This function is definitely steeper compared to the Sip ('!')- 
distributions which have been derived in Fendt (1997b) and are 
used in the present paper Here, we see the limitation of our 
ansatz. A steeper profile for rotation law is not yet possible to 
treat with our code due to the lack of numerical resolution. The 
non-linear character of the GS equation becomes more prob- 
lematic due to the gradients in the ilp -source term. 

In summary, only a model including differential rotation 
ilF(^') may provide a connection between the asymptotic jet, 
the disk magnetic flux distribution and also the size of the cen- 




0.10 0.12 0.14 0.16 0.18 0.20 

X 

Fig. 7. Magnetic flux distribution 'I'(x) along the disk surface 
as determined from the asymptotic jet properties and the Kep- 
lerian rotation of the disk. 



tral object. With our model we have presented a reasonable first 
solution for a self-consistent treatment. 



4.4. Application to the M87 jet 

The jet of M87 shows superluminal motion clearly indicating a 
highly relativistic jet velocity (Biretta et al. 1999). Recent radio 
observations have been able to resolve the innermost region 
of the M87 jet formation region with 0.33 x 0.12 mas beam 
resolution (Junor et al. 1999), corresponding to 2.5 — 7.0 x 
lO^^cm. Assuming a central supermassive black hole of 3 x 
IO^Mq (Ford et al. 1994), this is equivalent to about 30i?s! 
The derived jet full opening angle is 60° up to a distance of 
0.04 pc from the source with a "strong collimation" occurring 
afterwards (Junor et al. 1999). 

We now apply our two-dimensional jet model to these ob- 
servations and compare the geometrical scales. Such a com- 
parison is not possible for e.g. self-similar models. From the 
observed radio profile resolving the inner M87 jet (see Fig. 1 
in Junor et al. 1999), we deduce a jet radius of about 120 
Schwarzschild radii. With this, the first important conclusion is 
that the ratio of jet radius to light cylinder radius must be def- 
initely less than the value of 100 which is usually assumed in 
the literature. A number value of 3-10 seems to be much more 
likely. Numerical models of two-dimensional general relativis- 
tic magnetic jets fitting in this picture were calculated by Fendt 
(1997a). These solutions, however, do not take into account the 
differential rotation f2F(^'). 

Junor et al. (1999) claim that the M87 jet radius in the re- 
gion "where the jet is first formed cannot be larger than" their 
resolution of 30 Rs- Our conclusion is that the expansion rate 
is limited in both directions. The new radio observations give 
a minimum value of 3. Theoretical arguments limit the expan- 
sion rate to the value of about 20, since the jet mass flow must 
originate outside the marginally stable orbit which is located at 
3 — 6i?s. Clearly, if the jet radius is really as small as observed 
in M87, general relativistic effects may vary the field structure 
in the jet formation region. 
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From our model solutions, we derive a light cylinder radius 
of the M87 jet of about 50 Rs- The value derived from Eq. (14) 
differs from that by a factor of two, but is biased by the un- 
known size of the disk radius a;disk- This parameter, however, 
does not affect the global solution. Considering the standard 
relativistic MHD theory, nothing special is happening at the 
light cylinder For a highly magnetized plasma wind the light 
surface corresponds to the usual Alfven surface which itself 
does not affect the flow of matter. Hence, the light cylinder is 
un-observable. 

Also the opening angle in our numerical solution is larger 
than the observed value by a factor of two. This cannot be due 
to projection effects since any inclination between jet axis and 
the line of sight will increase the observed opening angle. We 
hypothesize that a numerical model with a steeper profile for 
the iso-rotation parameter will give a smaller jet opening angle 
comparable to the observed data. This is not surprising, since 
the jet footpoint anchored in a Keplerian disk rotates faster than 
in our model. Nevertheless, comparing the collimation distance 
observed in the M87 jet and assuming a similar ratio of jet ra- 
dius to light cylinder radius as in our model with h = 0.5, we 
find good agreement. The collimation distance is 2i?o- 

In summary, we conclude that the example of the M87 jet 
gives clear indication that the light cylinder of AGN jets might 
not be as large as previously thought. Although our model does 
not fit the observed geometrical properties of the inner M87 jet 
perfectly, we find in general a close compatibility. 

5. Conclusions 

We have investigated the two-dimensional magnetic field dis- 
tribution in collimating, relativistic jets. The structure of the ax- 
isymmetric magnetic flux surfaces is calculated by solving the 
relativistic force-free Grad-Shafranov equation numerically. In 
relativistic MHD, electric fields become important in difference 
to Newtonian MHD. The simplifying assumption of the force- 
free limit has been applied as relativistic jets must be highly 
magnetized. 

The central point of our paper is the consideration of dif- 
ferential rotation of the foot points of the field lines, i.e. a 
variation of the iso-rotation parameter fipC^)- The underlying 
model is that of a magnetic jet anchored in an accretion disk. 
Two main problems had to be solved in order to calculate a 
two-dimensional field distribution: a) to determine the a pri- 
ori unknown location of the light surface, b) the proper treat- 
ment of the regularity condition along that light surface. The 
light surface is the force-free equivalent of the Alfven surface 
and provides a singularity in the Grad-Shafranov equation. We 
summarize our results as follows. 

(1) We find numerical solutions for the two-dimensional 
magnetic flux distribution connecting the asymptotic cylindri- 
cal jet with a differentially rotating disk. In our example solu- 
tions the asymptotic jet radius is about 2.5 times the asymptotic 
light cylinder radii. This is the first truly two-dimensional rel- 
ativistic solution for a jet magnetosphere including differential 
rotation of the iso-rotation parameter Of('I'). The physical so- 



lution, being characterized by a smooth transition across the 
light surface, is unique for a certain parameter choice for the 
rotation law Op. 

(2) The half opening angle of the numerical jet solution is 
about 60 degrees. Cylindrical collimation is achieved already 
after a distance of 1-2 asymptotic light cylinder radii along the 
jet axis. Differential rotation decreases the jet opening angle, 
but increases the distance from the jet origin where collimation 
is achieved. The "jet expansion rate", defined as the ratio of the 
asymptotic jet radius to the jet radius at the jet origin, is about 
10. 

(3) From the analytical treatment of the regularity condition 
along the asymptotic branches of the light surface we derive a 
general estimate for the jet opening angle. We find that the jet 
half opening angle is larger than 45° and increases for a steeper 
profile of the differential rotation Hp. 

(4) Our two-dimensional ansatz, in combination with the 
treatment of differential rotation, allows for a connection of 
the asymptotic jet solution with the accretion disk. Certain disk 
properties can be deduced from the asymptotic jet parameters. 
Examples are the disk toroidal magnetic field distribution, with 
a maximum at half of the disk radius and the angular momen- 
tum flux per unit time and unit radius. This is interesting as 
a boundary condition for accretion disk models. We find that 
most of the angular momentum is lost in the outer part of the 
disk. 

(5) Application of our model to the M87 jet gives good 
agreement quaUtatively. From our numerical solution we de- 
rive an asymptotic light cylinder of the M87 jet of about 50 
Schwarzschild radii. Collimation of the jet would be achieved 
after a distance of two asymptotic light cylinder radii from the 
source. This value is comparable with the observations, how- 
ever, the opening angle in our model is larger by a factor of 
two. 

In summary, we have presented the first global two- 
dimensional solutions for a relativistic jet magnetosphere tak- 
ing into account differential rotation of the jet footpoints. From 
our jet model we may determine certain physical quantities in 
the disk that are not possible to observe, as e.g. the angular mo- 
mentum flux distribution at the disk-jet interface. Comparison 
with the M87 jet shows that our model seems to be consis- 
tent with the observations, therefore allowing for a derivation 
of the collimation distance, the light cylinder radius and the 
jet expansion rate for that example. Clearly, such features as 
the time-dependent ejection of knots and the interaction pro- 
cess between disk, jet and central source cannot be answered by 
our approach. Time-dependent relativistic MHD simulations of 
the whole "star"-disk-jet system would be necessary, however, 
such codes are not yet fully developed. 
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Appendix A: numerical methods 

For the solution of the two-dimensional GS equation we ap- 
ply the method of finite elements as developed by Camenzind 
(1987) and Fendt et al. (1995). Differential rotation f^pl^-) im- 
plies two major complications for the numerical computation. 
The first one is the fact that position and shape of the light 
surface 13 = is not known a priori. Along the light surface 
the boundary condition is the regularity condition, which, how- 
ever, itself depends on the two-dimensional solution z). 
The second problem is the GS source term for the differential 
rotation, containing the gradient of the magnetic flux, |V^p. 
Compared to the case of rigid rotation, this introduces another 
(and stronger) non-linearity in the GS equation. Therefore, a 
fragile numerical convergence process can be expected. 

An additional complication is that our grid of finite ele- 
ments of second order may be inadequate for a calculation of 
monotonous gradients between the elements if the numerical 
resolution is too low. However, for appropriate numerical pa- 
rameters as grid size, element size and iteration step size, we 
were finally able to overcome these difficulties. 

A.l. Determination of the light surface 

Here we discuss the iteration procedure we use to determine the 
location of the light surface. Because the rotation law riF(^) is 
prescribed, the radius where the light surface, D = {), intersects 
the jet boundary, = 1, is known, 

xlI^- = 1) = i/ripl^- = 1). (A.l) 

However, the corresponding position in z-direction is not 
known. Some estimates can be made about shape and inclina- 
tion of the light surface in the limit of large radii (see Sect. 2.4), 
but a general solution is not yet known. 

We start the iteration procedure calculating the inner solu- 
tion (defined as the field distribution inside the light surface) 
with an outer grid boundary at x = 1 (for comparison see 
Fig. 1). This choice is equivalent to the light cylinder in the 
case of rigid rotation. For differential rotation the radius x = 1 
is defined as asymptotic light cylinder (for large z). For low z- 
values the boundary x = 1 is located inside the light surface 
= 1/^f(^)- Along this outer boundary (of the inner 
solution), we apply a homogeneous Neumann boundary con- 
dition. Usually, this implies that the field lines will cross that 
boundary perpendicularly. However, in our case the homoge- 
neous Neumann boundary condition transforms into the regu- 
larity condition if the boundary becomes equivalent to the sin- 
gular light surface. As shown in Fendt et al. (1995), this trans- 
formation applies "automatically" in our finite element code. 
This is due to the facts that (i) finite element code solves the 
integrated GS equation and (ii) the boundary integral, which is 
proportional to D = 1 — a;^fi|, vanishes along the light surface. 

With the GS solution of the first iteration step we estimate 
the deviation of the chosen outer boundary from the true light 
surface by calculating D ~ 1 — x^f^K^I'). For the lowest z- 
value prescribed, we know that D ~ 1 — x^f2|(^ = 1). Then, 



the outer grid boundary (x, z) is slowly moved to a larger ra- 
dius with Ax ~ D{x, z)^. As a consequence of the different 
numerical grid, the field distribution will change. The value of 
D will, however, decrease. This procedure is repeated until D 
is below a certain limit, D ~ 0. Having obtained the solution 
inside the light surface, that field distribution is taken as inner 
boundary condition for the outer solution. 
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